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The purpose of a phase-preserving linear amplifier is to make a small signal larger, re- 
gardless of its phase, so that it can be perceived by instruments incapable of resolving the 
original signal, while sacrificing as little as possible in signal-to-noise. Quantum mechanics 
limits how well this can be done: a high-gain linear amplifier must degrade the signal-to- 
noise; the noise added by the amplifier, when referred to the input, must be at least half 
a quantum at the operating frequency. This well-known quantum limit only constrains the 
second moments of the added noise. Here we derive the quantum constraints on the entire 
distribution of added noise: we show that any phase-preserving linear amplifier is equivalent 
to a parametric amplifier with a physical state for the ancillary mode; the noise added to 
the amplified field mode is distributed according to the Wigner function of the ancilla state. 



I. INTRODUCTION 

The study of quantum limits on linear amplifiers became important in the 1960s with the in- 
vention and use of masers as microwave amplifiers. Initial investigations [IMS] led to the realization 
that quantum mechanics requires all phase-preserving linear amplifiers to add noise, thereby de- 
grading the signal-to-noise ratio of the input signal. For a high-gain linear amplifier, the minimum 
amount of added noise, when referred to the input of the amplifier, is equivalent to half a quantum 
at the operating frequency [2l[6]. For a quantum-limited input signal, this means a doubling of the 
input signal's zero-point noise and a halving of the input signal-to-noise ratio. 

This fundamental quantum limit is expressed formally as a bound on the second moment of the 
noise added by a phase-preserving linear amplifier. A comprehensive review article reprises the 
development and elaboration of this fundamental quantum limitation on the operation of linear 
amplifiers [7j. In recent years, microwave-frequency amplifiers, based on the Josephson effect, 
have very closely approached the fundamental quantum limit on second-moment added noise [8l- 
110] . In the meantime, workers in quantum optics have formulated techniques for determining 
photon correlation functions without using photon counting, instead using the linear detection 
that at optical frequencies comes from homodyne detection |llj . Researchers working with linear 
amplifiers at microwave frequencies have refined and elaborated these techniques into methods for 
determining the noise properties of signals input to a linear amplifier and of the added amplifier 
noise |12H14j . These methods have been used to determine moments of amplifier noise well beyond 
second moments \13\ [T5] , to measure photon correlation functions of input microwave signals |16| , to 
do quantum tomography on itinerant (wave-packet) microwave photons |17j . and to study squeezing 
of microwave fields |18| [T9] . 

All these developments motivate an investigation of quantum limits on all moments of the added 
noise or, equivalently, on the entire distribution of added noise. Second moments are sufficient to 
characterize the added noise if it is Gaussian; measuring higher moments allows one both to check 
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the Gaussianity of the added noise and to characterize the performance of linear amplifiers more 
thoroughly. Here we consider the case of phase-preserving amplification of a single bosonic mode, 
which we call the primary mode and which has annihilation operator a. We characterize the input 
and output noise in terms of symmetrically ordered moments of a and or, equivalently, in terms 
of symmetrically ordered moments of the input and output quadrature components. With this 
convention, the noise is described completely by the input and output Wigner functions of the 
primary mode. 

We show here that regardless of how a phase-preserving linear amplifier is realized physically, 
it is equivalent to a parametric amplifier, i.e., an amplifier in which the primary mode undergoes 
a two-mode squeezing interaction with a single ancillary mode, which has annihilation operator b. 
The strength of the parametric interaction determines the amplifier's gain, and the noise added 
by the amplifier is distributed according to the Wigner function of the ancillary mode's initial 
state a. Characterizing completely the noise properties of a phase-preserving linear amplifier thus 
amounts to giving the initial state of this effective ancillary mode, even though the amplifier might 
be nothing like a parametric amplifier. A quantum-limited (ideal) linear amplifier corresponds to 
the case where a is the vacuum state. 

We begin by reviewing in Sec. II A the simple input-output relation that leads to the second- 
moment constraint on added amplifier noise. An ideal linear amplifier saturates the second-moment 
constraint and has Gaussian noise. In Sec. II B we give a stick-figure pictorial representation of the 



input and output noise in terms of contours of the popular quasidistributions for a field mode [20- 
[23] . the Glauber-Sudarshan P function [23H26] . the Wigner function W [27], and the Husimi Q 
distribution [28 1, and in Sec. II C we consider four generic models for an ideal linear amplifier. In 
Sec. Ill we develop a general mathematical description of a linear amplifier that adds arbitrary 
noise. The amplifier is described in terms of a linear map that takes the input state to the output 
state; this amplifier map must be completely positive \29\ [30] to correspond to a physical linear 
amplifier. Section [IV] formulates and proves our main result: the requirement of complete positivity 
implies that any linear amplifier is equivalent to a parametric amplifier with a physical initial state 
a for the ancillary mode. Section [V] considers examples of nonideal amplifiers, including unphysical 
ones, and Sec. VI uses our main result to spell out the quantum limits on higher moments of the 
added noise. Section VII sums up and briefly sketches future work. The manipulations necessary 
to relate various kinds of moments are given in an Appendix. 



II. QUANTUM-LIMITED PHASE-PRESERVING LINEAR AMPLIFIERS 
A. Quantum limit on second moment of added noise 

The setting for our investigation is a single bosonic mode, called the primary mode, which is to 
undergo phase-preserving linear amplification. The primary mode has annihilation and creation 
operators, 

a = -j=(xi + ix 2 ) , (2.1) 

^ = —(x 1 -i X2 ). (2.2) 

in these expressions, the rapid oscillation at the modal frequency cj has been removed, and x\ and 
X2 are the Hermitian quadrature components of the mode. The creation and annihilation operators 
obey the canonical commutation relation, [a, a*] = 1 (equivalently, [xi,X2] = i)- This implies an 
uncertainty principle, (Ax 2 ) (Ax 2 .) > 1/4, where A denotes the difference between an operator and 
its expectation value, Ax = x — (x), and hence (Ax 2 ) is the variance of x. 
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We can think of the signal as being carried by a single-mode field 

E(t) = \(ae~ iult + a j e~ iujt ) = —= (an cos wt + x 2 sin ut) . (2.3) 
2 v 2 

The annihilation operator is a complex-amplitude operator for the field; the expectation value of 
the field, (E(t)) = Re((a) e~ lujt ) , oscillates with the amplitude and phase of (a). The variance of i£ 
characterizes the noise in the signal; for phase-insensitive noise, this variance is constant and given 
by (AE 2 ) = l(\Aa\ 2 ), where 

(|Aa| 2 ) = ^<AaAa f + Aa+Aa) = (|a| 2 ) - |(a)| 2 (2.4) 

is the symmetric variance of a. Here we use the notation | a \ 2 = \{aa) + a) a) for the symmetric 
product of a and a) [31] . 



The symmetric variance (2.4) obeys an uncertainty principle, 



(|Aa| 2 ) = \{(Ax 2 ) + (Ax 2 )) > (Ax 2 )" 2 (Ax 2 )" 2 > \ . (2.5) 

The lower bound is the half-quantum of zero-point (or vacuum) noise. The first inequality is 
saturated if and only if the noise is phase insensitive, i.e., (|Aa| 2 ) = (Ax 2 ) = (Ax 2 ), the second if 
and only if the quadrature uncertainties have minimum uncertainty product. Both inequalities are 
saturated if and only if the mode is in a coherent state \a) = D(a, a)|0), where 

D(a, a) = e aa ^ a * a = e i{a ^~ a ^ , a = ^=(a x + ia 2 ) , (2.6) 

v2 

is the displacement operator for mode a. We use a two-slot notation for the displacement operator, 
partly so as to identify the mode the displacement operator pertains to and partly so that by 
putting a c-number in both slots, as in D(/3,a) = e a @*~ a *P = e*( a2 ^ 1_ai ^ 2 ), we have a convenient 
notation for two-dimensional Fourier transforms |31| . 

The objective of phase-preserving linear amplification is to increase the size of an input signal 
by a (real) multiplicative amplitude gain g, regardless of the input phase, while introducing as little 
noise as possible. The amplification of the input signal can be expressed as the transformation 

(flout) = g{ata) (2.7) 

of the expected complex amplitude. A perfect linear amplifier would perform this feat while 
preserving the signal-to-noise; in the Heisenberg picture, the primary mode's annihilation operator, 
not just its expectation value, would transform from input to output as 

flout = gain ■ (2.8) 

The second-moment noise would be amplified by the power gain g 2 , i.e., (|Ao out | 2 ) = g 2 (\Aai n \ 2 ) . 
The amplifier's output would be contaminated by the same noise as the input, blown up by a factor 
of g 2 , but the amplification process would not add any noise to the amplified input noise. 

Alas, quantum mechanics prohibits free lunches: there are no perfect phase-preserving linear 



amplifiers; the transformation (2.8) does not preserve the canonical commutation relation and thus 
violates unitarity. Physically, this is the statement that amplification of the primary mode requires 
it to be coupled to other physical systems, not least to provide the energy needed for amplification; 
these other systems, which can thought of as the amplifier's internal degrees of freedom, necessarily 
add noise to the output. This physical requirement is expressed in an input-output relation [31 16], 

a ou t = gain + L ] , (2.9) 
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where the added-noise operator L is a property of the internal degrees of freedom. One usually 
assumes that (Lt) = so as to retain the expectation-value transformation (|2.7|). Preserving the 
canonical commutation relation between input and output requires that 



[L, Lt] 
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(2.10) 



which implies an uncertainty principle, 



|AL| 2 )> l -(g 2 -l 



(2.11) 
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FIG. 1: Input signal E- m (t) = cos cot (blue) and output signal E out (t) = g cos Lot (red) for a single-mode, 
phase-preserving linear amplifier with amplitude gain g = 4. Input and output have the same phase. The 
input complex amplitude, (oi n ) = (xi)/y/2 — 1, is amplified to become the output complex amplitude, 
(flout) = g- The (phase-insensitive) noise on the signal is represented by smearing out the mean signal into 
a band with vertical height equal to the uncertainty in the field, (AE 2 ) 1 / 2 = dAa] 2 ) 1 / 2 /s/2; for the ideal 
linear amplification with coherent-state input shown here, the height of the input band is |, and the height 



of the output band is y (g 2 — |)/2. The phase-space diagrams depict the same input and output as do the 
temporal plots. The mean complex amplitude is represented by an arrow, and the noise by a circle whose 
diameter is equal to the height of the band in the temporal plot. The temporal plots can be obtained by 
rotating the phase-space stick and ball about the origin and projecting onto the real axis. 

The amplifier should be prepared to receive any input in the primary mode, without having 
any idea what that input is going to be. This places the restriction that the primary mode and the 
internal degrees of freedom cannot be correlated before amplification. The total output noise is 
then the sum of the amplified input noise and the noise added by the internal degrees of freedom: 

(|Aa out | 2 )= 5 2 (|Aa in | 2 ) + (|AL| 2 ). (2.12) 



The added noise is constrained by the uncertainty principle (2.11), which together with Eq. ( |2.5| ), 
places a lower bound on the output noise: 

(|Aa out | 2 ) > g 2 - \ . (2.13) 

There are, of course, states for which the output noise is smaller — indeed, as small as half a 
quantum — but these require that the primary mode and the amplifier's internal degrees of freedom 
be correlated at the input. Figure [T] illustrates the amplification of the field E and introduces the 
traditional ball-and-stick phase-space diagrams that are used to depict the second-moment noise. 
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Referred to the input, the output noise takes the form 

= <|W) + = (|A ,„| 2 > + A ■ (2.14) 

9 9 

The second- moment added noise, referred to the input, is called the added-noise number A [6]. 
It provides a natural, dimensionless characterization of an amplifier's performance, and it is con- 
strained by quantum mechanics to satisfy 

g 2 ~ 2 V g 2 J Q^oo 2 V ' 

Amplifier performance is often characterized by noise temperature T n or noise figure F. The 
noise temperature is defined as the temperature required to account for all of the output noise 



referred to the input, i.e., to account for the noise (2.14); when the input is quantum limited, i.e., 
a coherent state with (|Aai n | 2 ) = i, the noise temperature is related to the added-noise number 
by A = (e riu >/ k B T n _ The noise figure is the ratio of the input signal-to-noise ratio to the 

output signal-to-noise ratio; for quantum-limited input, the noise figure is given by F = 1 + 2A. 
For quantum- limited input, the noise temperature and noise figure satisfy 

ksT n 1 1 1 

~lv7 = ln(l + „4- 1 ) - ln[(3-l/<? 2 )/(l-l/<? 2 )] ^ hi3 ' ( ' ' 

F = 1 + 2A > 2 - 1/g 2 — >2; (2.17) 

an amplifier that operates far from the quantum limit has A = F/2 = kBT n /hoj 3> 1. In the limit 
of high gain, a phase-preserving linear amplifier adds at least half a quantum of noise to the input 
noise, and as a consequence, for a quantum-limited input, the signal-to- noise ratio is degraded by 
at least a factor of two. 

All three of these measures of amplifier performance are afflicted by the residual gain dependence 
1 — 1/g 2 . The noise temperature and noise figure have the additional annoyance that, as they depend 
on the input noise, they are not solely properties of the amplifier. Finally, the noise temperature 
is not even linear in the added noise for amplifiers operating near the quantum limit. We prefer in 
this paper to deal with an added-noise number that has all the gain dependence removed: 

(|AL| 2 ) A 1 , \ 

At = M „ 17 = > - . 2.18 

g 2 - 1 1 - 1/g 2 " 2 { ' 

The second-moment quantum limit is A\ > 1/2. The subscript indicates that Ai is the first in a 
sequence of added-noise numbers. We introduce added-noise numbers for all moments of the added 



noise in Sec. VI and consider the limits imposed by quantum mechanics on moments of all orders. 



B. Ideal linear amplifier 



An ideal linear amplifier saturates the second-moment bound ( 2.18 ) . The added noise in this case 



is necessarily Gaussian, as we show in Sec. IV; there are no constraints on an ideal linear amplifier 
beyond this second-moment bound. It is instructive to introduce here a pictorial representation of 
the amplified input noise and the added noise for the case of an ideal linear amplifier acting on a 
quantum-limited (coherent-state) input. This allows us to discuss the several perspectives provided 
by the various ways of ordering creation and annihilation operators. Since the noise is Gaussian, 
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Q distribution 



FIG. 2: Ball-and-stick phase-space depictions of input and output states for an ideal linear amplifier whose 
input is a coherent state |/3) with \j3\ — 1 and which has amplitude gain g — 4, giving the output state 
a mean that lies on a circle of radius g\a\ — 4 centered at the origin. The three cases correspond to 
normal ordering (Glauber-Sudarshan P function), symmetric ordering (Wigner W function), and antinormal 
ordering (Husimi Q distribution) . The input (blue) and output (red) noise are characterized by the variance 
E 2 of a: the radius of the noise circle is chosen to be E/2-\/2. For normal ordering, Ep = (Aa^Aa); for 
symmetric ordering, Y? w = (|Aa| 2 ) = \{AaAa> + Aa^Aa) (E w /2\/2 = ^(AE 2 ) 1 / 2 ); and for antinormal 
ordering, Eg = (AaAa 1 '). At the input, Ep = 0, Efy = \, and Eg = 1; at the output, Ep = g 2 - 1, 
E 2 ^ = g 2 — ~, and Eg = g 2 . For each ordering, blue is used for the input state and for the amplification of 
the input mean and noise by a factor of g; red depicts the total noise at the output. For symmetric ordering, 
the output noise consists of the amplified input noise, with variance \g 2 1 and added noise, with var ianc e 



\{g 2 — 1), which degrades the output signal-to- noise ratio relative to the perfect amplification of Eq. (2. 



For normal ordering, the input is represented by a (noiseless) point, which is amplified to another point; 
all of the output noise, which has variance g 2 — 1, appears to be added noise. The flip side is antinormal 
ordering, where the input has an extra half-quantum of noise relative to symmetric ordering; the amplified 
input noise, which has variance g 2 , accounts for all of the output noise, and there appears to be no added 
noise. Even for the modest gain used here, it is difficult to distinguish the sizes of the output noise circles 
for the three ordcrings, because the output noise dwarfs the half-quantum difference between the orderings. 



the pictorial representation can be simplified to ball-and-stick figures, like those in Fig. [TJ which 
depict only the first and second moments. 

A Gaussian phase-space distribution that has phase-insensitive noise has the form 
e -\a-$\ /£ /■kT? , where [3 is the mean value of a and T? is the variance of a; this variance can be 
calculated using several orderings, which go with different quasidistributions for the field mode [20f - 
123] . Up till now, we have used the symmetrically ordered variance = (|Aa| 2 ) = 2(AE 2 ), which 
goes with the symmetric ordering of the Wigner function [27J and which we thus denote with a sub- 
script W. If we use normal ordering of a and a\ the appropriate quasidistribution is the Glauber- 
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Sudarshan P function [2~H - [2"6] . and the normally-ordered variance is T? p = (AatAa) = — |. 
Likewise, if we use antinormal ordering, the appropriate quasidistribution is the Husimi Q distri- 
bution [28], and the antinormally-ordered variance is £q = (AaAa^) = + |. 

For an ideal linear amplifier with coherent-state input, we depict the input and output in Fig. [2j 
for normal, symmetric, and antinormal ordering of the noise variance. We use circles centered at 
the mean value, with the radius S/2\/2 representing the size of the noise; the stick of Fig. [I] is 
omitted as redundant. The multiple of £ we use for the radius of the noise circle is chosen so that 
the noise circles fit within the figure without overlapping. One other ingredient appears in Fig. [2j 
the amplified input noise is obtained by expanding the input-noise circle by a factor of g. For the 
case of symmetric ordering, this amplified input noise shows what the output noise would be for a 
perfect linear amplifier, as in Eq. (2.8). The other two orderings give different perspectives on the 
amplifier noise, discussed in the figure caption. 

Of particular interest is the Husimi Q distribution, which for a modal state p is given by 
Q p {a) = (a\p\a)/ir. The stick-figure depiction of Fig. [2] shows that from the perspective of the 
antinormally-ordered variance, all the output noise in an ideal linear amplifier is amplified input 
noise, with no added noise at all. Indeed, Fig. [2] suggests that the output Q distribution of an ideal 
linear amplifier is a scaled version of the input Q, scaled by the gain g, i.e., Qout(P) = Qm(P/g)/g 2 , 
and this turns out to be true for arbitrary input states, as we show shortly. There is a reason why 
there is no added noise in the Q-distribution picture: the nonnegative Q distribution describes the 
statistics of a quantum-limited, simultaneous measurement of both quadrature components, x\ and 
X2] the scaling of the Q distribution from input to output says simply that relative to quantum- 
limited measurements of the signal at both the input and the output, there is no degradation of 
signal-to-noise ratio between the input and the output of an ideal linear amplifier. We stress this 
conclusion: relative to the best measurements one can make, there is no loss of signal-to-noise in 
amplifying a signal. The reduction in signal-to-noise occurs in symmetrically ordered moments, not 
in the antinormally ordered moments that apply to simultaneous measurements of both quadrature 
components. 

The goal of this paper is to go beyond the Gaussian noise of an ideal linear amplifier. Thus we 
will need to move beyond the ball-and-stick diagrams of Fig. [2] and plot the entire distributions of 
input and output noise. In doing so, we will make use of all three operator orderings, moving freely 
among them as we find one or the other better serves our purpose. We find it most convenient to 
formulate our general mathematical description of a phase-preserving linear amplifier in terms of 
the P-function picture, though we quickly generalize to all three orderings. Before turning to that 
task, which occupies Sec. Ill, we review several models of an ideal linear amplifier. 

We also note complementary work on the quantum noise limits for operational amplifiers 
[32j 133] . An operational amplifier takes in the voltage, say, at the end of a transmission line and 
outputs a voltage into another transmission line. Since the input and output voltages both consist 
of incident and reflected waves, an op-amp does not have the simple input-output structure that 
is the basis of our work here and most previous work on linear amplifiers. 



C. Models for ideal linear amplifiers 

In this subsection we consider models for ideal linear amplifiers. While detailed microscopic 
models of particular amplifiers and their noise sources (see, e.g., [Ml E5]) are important, they 
are not central to our analysis. Here we survey four generic models of a linear amplifier, to 
build intuition about the fundamental physical processes that account for amplifier noise and to 
provide context for our subsequent study of general quantum constraints on the performance of 
linear amplifiers. 



s 



1. Parametric amplifier 



The simplest model of an ideal linear amplifier is provided by a parametric amplifier |36H38| . 
The primary mode a interacts with an ancillary mode b = (y\ + iy2)/V%, which is initially in the 
vacuum state. The total Hamiltonian, 

H = hw{a)a + tfb) + ihK(abe 2iujt - a^b^e~ 2iujt ) , (2.19) 

has an interaction term that describes pairwise creation or destruction of quanta in the two modes. 
This pairwise creation or destruction is accompanied by destruction or creation of a quantum in a 
pump mode that has frequency 2u. The pump mode does not appear in the Hamiltonian because 
it is excited into a high-amplitude coherent state and thus is essentially classical. Its amplitude 
contributes to the coupling strength k, and its time dependence gives the e ±2 * w * explicit time 
dependences in the Hamiltonian. 

Transforming to the interaction picture that removes the free Hamiltonians of the two modes, 
the interaction part of the Hamiltonian assumes the form 

il/ = ihn(ab - attf) , (2.20) 

which can be integrated to give an evolution operator 

C7>(*) = e~ iHlt/h = exp[r(a6 - a)b ] )] = exp[ir( Xl y 2 + x 2 yi)] = S(r) , r = nt , (2.21) 

where S(r) is the two-mode squeeze operator [39-41J. In the Heisenberg picture, the primary mode's 
annihilation operator undergoes the transformation 



flout 



S^aS = a cosh r — b^ sinh r = ga — b^ \/ g 2 — 1 , (2.22) 



i.e., the amplifier input-output relation (2.9) with gain g = coshr and noise operator L = 
—b sinh r = —b\J g 2 — 1. If the ancillary mode begins in the vacuum state |0), the noise opera- 
tor saturates the second-moment bound (2.18), and we have an ideal linear amplifier. 

We can also describe the evolution in the interaction picture. If the primary mode has initial 
state p and the ancillary mode begins in vacuum, the state of the primary mode after amplification is 

£{p) = tv b [S{r) P ® |0><0|5^(r)) . (2.23) 

Here E is the trace-preserving quantum operation (a completely positive map) that describes how 
the state of primary mode transforms from the input to the output of the amplifier. Figure [3^a) 
gives the simple quantum circuit for an ideal parametric amplifier. 
The factored expression for the squeeze operator [4*01 
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can be used to eliminate the ancillary mode from the ideal-amplifier quantum operation ( |2.23 ). 
In a first approach, we find the partial matrix element of S(r) between a coherent state |/3) and 
vacuum for the ancillary mode: 



{P\S(r)\0) = e — g-^ a e-V9^P*^ = . (2.25) 



Taking the trace in Eq. (2.23) in the coherent-state basis of the ancillary mode, 

l2f. 



; (p ) = J*l (f}\S\0)p(Q\Sfi\[}) = J d 2 pA pP A\ (2.26) 
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(a) 



P 



|0> 



S(r) 



£(p) = tr b (S(r)p®\0)(0\S^r)) 
Joint output state S(r)p® |0}(0|S%) 



(b) 



S(r) 



£(p) = tr b (S{r)p®a&(r)) 
Joint output state S(r)p (g) aS^(i 



FIG. 3: (a) Ideal phase-preserving linear amplifier realized as a parametric amplifier. The primary mode a 
interacts with an ancillary mode b, which begins in the vacuum state |0), via a two- mode squeezing interac- 
tion. The output state £{p) of the primary mode is amplified with an amplitude gain g = coshr and with 
the minimum amount of added noise permitted by quantum mechanics, (b) Parametric amplifier with an 
arbitrary initial state a for the ancillary mode. The main result of this paper is that any phase-preserving 
linear amplifier, no matter its physical realization, is equivalent to a parametric amplifier with some physical 
initial state a for the ancillary mode. 



gives a Kraus decomposition of £; the operators Ap are called Kraus 
operators |29} 130] . One can use this to find the output Q distribution for an ideal linear amplifier. 
As promised above, the output Q distribution is a scaled version of the input Q distribution: 

Q out ( a ) = -(a\£(p)\a) = . (2.27) 

A different Kraus decomposition of £ [32] follows from partial matrix elements of the squeeze 
operator in the number basis of mode b: 

(— Tin ( n 2 _ i \n/2 

<ra|S(r)|0> = 1 ' [9 ' ff-° ta (at)" = A n . (2.28) 
9 Vnl 



Taking the trace in Eq. (2.23) in the number basis of the ancillary mode gives the Kraus decom- 
position 

oo oo 

£{p) = £>|5|0>p<0|5t|n) = A nPA . (2.29) 

n=0 n=0 



2. Inverted-oscillator model 



Another model for an ideal linear amplifier, due to Glauber [43J, uses a primary mode a and 
an ancillary mode b. The two modes have the same frequency u, but the ancillary mode is an 
inverted oscillator (sometimes called a negative- mass oscillator [44H46J ) . The inverted oscillator 
has an upside-down Hamiltonian, —hojb^b. Since its energy levels run down instead of up, the 
inverted oscillator is a source of energy; when a quantum is created in the inverted oscillator, the 
oscillator emits energy hu. The Hamiltonian of the two modes is 



H = huj(a ] a - b ] b) + ihn{ab - a ] b ] ) 



(2.30) 



10 



Transforming to an interaction picture that removes the free Hamiltonians gives the interaction 
Hamiltonian (2.20) of a parametric amplifier; the subsequent discussion of linear amplification is 
thus identical to that for a parametric amplifier. Indeed, the only difference between a parametric 
amplifier and the inverted-oscillator model is that the inverted oscillator does not need a pump 
to balance the energy books; creation of a quantum in the inverted oscillator provides the energy 
needed to create a quantum in the primary mode. 



3. Linear- amplifier master equation 



Our third model is an elaboration of either a parametric amplifier or the inverted-oscillator 
model. The single ancillary mode is replaced by a field, which is initially in the vacuum state. The 
instantaneous temporal field modes interact with the primary mode via a parametric interaction 



like Eq. (2.20); the result is the master equation for an ideal linear amplifier |47j . Models of this 
sort, based on coupling the primary mode to a sequence of inverted oscillators, have been developed 
by several authors 

This approach starts with an (interaction) Hamiltonian 

/•oo 

H = ih^j / dr (ab T - a ] b\)5(T - t) = ih^j{ab t - o)b\) . (2.31) 
J o 

The operators b T and bt are continuum annihilation and creation operators for the instantaneous 
field modes, obeying the canonical commutator [6 r ,fcL] = 5(t — r'). The parameter r labels the 
field modes and specifies the time t = r at which a field mode interacts with the primary mode. 
It is easy to derive the Heisenberg-picture equations of motion: 



da 

~dt 
db T 

~dJ 



1 

ih 
1 

ih 



[a,H] -- 
[br,H] 



-Vjb\(t) , 



The solution of Eq. (2.33), 



b T (t) = b T (0) - ^a\T)Q(t - t) , 



(2.32) 
(2.33) 

(2.34) 



where O(i) is the unit step function, with value \ at t = 0, can be plugged into Eq. (2.32) to give 



da 
~dt 



i 7 a(t)-v^6l(0) 



whose solution is 



Here 



a{t) = g{t)a{0) + L\t) . 



g(t) = 



is the amplitude gain that applies if the interaction is turned off at time t, and 



L{t) = ~Vl / dr g(t - r)6 T (0) 
J o 



(2.35) 



(2.36) 



(2.37) 



(2.38) 
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is the added-noise operator. It is easy to verify that L{t) and L'(t) satisfy the commutator (2.10), 
as required by unitarity. One can think of all the added noise as coming from a single, discrete, 



wave-packet mode, whose annihilation operator is b = L/y' g 2 — 1. 

The easiest way to derive the corresponding master equation is to discretize the field modes 
into wave packets, each of which lasts a short time At: 



'At 



drb T . 



(2.39) 



2-1 



Here tj = jAt. The primary m ode in teracts sequentially with these discretized modes, according 
to the interaction Hamiltonian (2.20) with coupling constant k = v/7/ At. The interaction of the 
primary mode with the jth discrete mode changes the state of the primary mode according to 



p(tj) = tr 6 . (5(KAt)p(^_i) ® \0){0\S\nAt)) 



(2.40) 



Expanding the squeeze operators to second order and rewriting in terms of the change in p through 
the jth interaction gives 



— — = -K 2 At{2a}pa — aa) p — paaJ) . 



At 2 

We now take the limit At — > and k — > 00, with n 2 At = 7 held constant, obtaining 



dp 



1 



(2a* 



pa — aw p — paa^) . 



(2.41) 



(2.42) 



This is the (ideal) linear- amplifier master equation |47j . 

Solving the master equation is easy. Using the standard rules |20j . translate Eq. (2.42) to a 
partial differential equation for the Q distribution: 







d 7 
8t + 2 



da da* 



Q(a, a*;t) . 



(2.43) 



The solution, 



Q(a, a*;t) 



a 



a 



9 2 {tr\g{ty g(t) 



;0 



(2.44) 



where g(t) is the amplitude gain (2.37), shows again that the input-output transformation for an 
ideal linear amplifier is a rescaling of the input Q distribution by the amplitude gain. 



4-. Measurement-based model of linear amplification 

The last model we consider explores the connection between ideal linear amplification and 
quantum-limited simultaneous measurements of both quadrature components [50]. The strategy 
for amplification is to measure both quadrature components of the primary mode, x\ and X2, and 
then to create an amplified coherent state \ga), where a = {a\ + ia^/v^ is determined by the 
measurement outcomes, a\ and «2- Discarding the outcomes introduces an average of the output 
coherent states over the probability distribution for the measurement outcomes. The result is a 
linear- amplification process, with the added noise due to the noise that accompanies a simultaneous 
measurement of x\ and %2- This is a silly strategy for making a linear amplifier, because the point 
of linear amplification is to make a small signal accessible without the need for quantum-limited 
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measurements. Nonetheless, this approach is instructive in highlighting the connection between 
quantum-limited amplification and quantum-limited measurements of the quadrature components. 

A quantum-limited simultaneous measurement of both quadrature components is described by 
coherent-state projectors. The Kraus operators for such a measurement are 



1 



7T 



a)(a| 



(2.45) 



The trace-decreasing quantum operation for outcome a is K a pKa = Q p (a)\a)(a\; i.e., the mea- 
surement statistics are given by the Q distribution, and the output state is the coherent state that 
corresponds to the measurement outcomes. The quantum-circuit diagrams in Fig. [4] summarize 
pictorially the results of the algebraic contortions in the following analysis. 



Our goal is to implement the Kraus operators (2.45) in an ancilla model. To do so, we introduce 
two ancillary modes, b = (yi + iy2)j\[2 and c = (z\ + iz2)/V2, which serve as meters that record 
the measurement results. For mode b, we introduce ^-normalized eigenstates of the quadrature 
components, Vj\Pj) - 
j = 1,2, with (71I72) 

in the first quadrature, y±, of mode 6; thus we identify a\ with f3\. Similarly, the result a 2 of 
measuring x 2 is recorded in the second quadrature, z 2 , of mode c; thus we identify a 2 with 72. 
We begin the analysis by using [201 [2H ED] 



j = 1,2, with = e^/v^; similarly, for c, z 3 \ lj 



= *7l 72 



/ \/2tt. In the ancilla model, the result a\ of measuring x\ is recorded 



|a)(a| 



d 2 fi 



7T 



d 2 n 



7T 



e- M2/2 D(a,fi)D{fi,a) 



to manipulate K a into the form |51j 

1 f dfcdji 



-(/3|+7^)/4 -i(/3 2 xi-7ix 2 ) „«(/3i/3 2 -7i72) 



2TT 



(2.46) 



(2.47) 



Here we relabel the measurement outcomes as described above, i.e., a\ = /3i and 02 = 72, so 
a = (Pi + i72)/v / 2, and we set \i = —(71 + i^/V^- We let the initial state \4>) of the ancillary 
modes be specified by the wave function (/32,7i|</>) = e~^ 2+7l ^ 4 /v / 27r. This initial state can be 
written as 



|$ = Si(r)®Sl(r)|0,0>, 
where |0, 0) is the vacuum state of modes b and c and 

Si(r) = exjp[lr(b 2 - 6 t2 )] = exp[^r(yiy 2 + 2/22/1)] 



(2.48) 



(2.49) 



is the single-mode squeeze operator for mode b (similarly for mode c) |3H [391441] . and where the 
squeeze parameter here corresponds to squeezing by a factor of 2, i.e., e 2r = 2. In the state \(p), 
the four quadrature components are uncorrelated and have variances 



(Ay?) 



(Ay 2 ) = (Az 



-2r 



1 

4 ' 



,2r 



1 



(2.50) 



We can now write the Kraus operator (2.47) in the form 



K a = V2 d(3 2 d 7l (/3i,72|/3 2 ,7i)(/3 2 ,7i|e- %2a;i - Zi:C2) |0) = v^i, 72^1! 



(2.51) 
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(a) |0) 



0) 



Si(r) 



S\(r) 



— e 



-i(y2Xi-zix 2 ) 




yi 




=*(S'-l)(72a;i-^iX2) 



72 



(b) |0) 



|0> 



Si(r) 



-i{y2Xi-ziX2) 



u 2 

J{g-l)(z2xi-yix2) 




01 



III 




72 



Z2 



FIG. 4: (a) Quantum circuit for measurement-based amplification of primary mode a = {x\ + ix 2 ) /v2. 
(b) Circuit for the equivalent coherent scheme. Filled circles represent control on the first quadrature 
of a mode; open circles represent control on the second quadrature. In both cases, the initial controlled 
operation is a sensing step in which the two quadratures, x\ and X2, of mode a are written onto ancillary 
modes b — (yi + iy 2 ) / V2 and c = [z\ + iz 2 )/v / 2; X\ is recorded in yi, and x 2 is recorded in z 2 . This sensing 
step is followed by a feedback operation from b and c to amplify the complex amplitude of a. In (a) the 
feedback operation is classical and based on the results, fi\ and 72, of measuring y\ and z 2 (measurements are 
represented by rounded boxes); discarding the measurement outcomes yields a linear-amplification process 
on the primary mode. In (b) the feedback operation is coherent; the final measurement is irrelevant to the 



amplification and can be omitted. The amplified output state (2.57) of mode a is obtained by tracing out 
modes b and c. The equivalence of the two circuits, demonstrated in the text, is obvious from the circuit 
diagrams as an example of the principle of deferred measurement [30] [55] ■ The single-mode squeezers applied 
initially to the ancillary modes prepare the input state \<p) of Eq. (2.48). When the squeeze parameter r 
is chosen to squeeze by a factor of two (e 2r = 2), the model is based on the Arthurs-Kelly simultaneous 
measurement of X\ and x 2 [SOJ, which is described by coherent-state projectors, but the result is not quite 
an ideal linear amplifier. When the squeeze parameter is specified by e 2r = 2[g — l)/(g + 1), the result is 
an ideal linear amplifier, based on a slightly different model of simultaneous measurement of x± and x 2 ; this 
realization of an ideal linear amplifier is different from a parametric amplifier. 



where 

Ui = exp[-i(y 2 x 1 - Zix 2 )] ■ (2.52) 

This demonstrates one way to implement K a : let the primary mode interact with the two ancillary 
modes via an instantaneous interaction H = h{y 2 x\ — zix 2 )8(t); then measure y\ and z 2 on the 
ancillary modes, getting results /3± and 72, giving a Kraus operator Kp ltl2 = K a /y2. The y/2 
comes from the change in integration measure in going from a = (/?i + i^ 2 )/y/^ to j3\ and 72. The 
interaction displaces the first 6-quadrature, yi, by x\ and the second c-quadrature, z 2 , by x 2 ; the 
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measurements of y\ and z 2 read out these displacements, contaminated by the uncertainties in y\ 
and z 2 , and leave the primary mode in the coherent state \a) corresponding to the measurement 
results, a\ = j3\ and a 2 = 72- 

To turn this into an amplifier model, we use the result of the measurement to displace the 
primary mode by (g — l)a. The entire procedure is then described by Kraus operators 

K' a = D(a, (g - l)a)K a = -^\ga)(a\ . (2.53) 

\/7T 



Writing this in terms of quadrature components of the ancillary modes gives 

K'^ 2 = K'JV2 = e^-V^-^X^MUA® = , (2.54) 

where 

U 2 = exp[i(g - l)(z 2 xi - 2/1X2)] • (2.55) 



In the first form in Eq. (2.54), there is a sensing interaction U\ and then measurements of y\ 
on mode b and z 2 on mode c; this is followed by an amplifying displacement D{a, {g — l)a) = 
e*(9- 1 )(72a;i-^i^2) Q f m ode a based on the measurement results fi\ and 72. Discarding the measure- 
ment outcomes leads to a linear amplifier. In the second form, there are two coherent interactions, 
first the sensing interaction U\ and then an amplifying feedback U 2 , and these two are followed 
by the measurements of yi and z 2 ; in this second form, discarding the measurement outcomes can 
be accomplished by omitting the closing measurement. These considerations are summarized in 
quantum-circuit diagrams in Fig. |4j 

The output state of this linear amplifier is 

£{p) = J d 2 aK' aP K£ = J d 2 a Qp ^J 9 \ a){a\ . (2.56) 

This can also be written as 

S(p) = tv biC (U 2 U lP 8) \(f>}(0\uM) . (2.57) 



It is easy to see from Eq. (2.56) that this is not quite an ideal linear amplifier: the output P 
function, not the output Q distribution, is a rescaled input Q distribution. This means that this 
amplifier adds two more units of vacuum noise than does an ideal linear amplifier. We can see 
this more directly — and see also how to convert to ideal linear amplification — by examining the 
input-output relation for the primary mode's annihilation operator, 

flout = u\ulaU 2 Ui =ga + ^g 2 - Id) . (2.58) 

Here we introduce a modal annihilation operator d = (si + is 2 )/^/2, where 



■''W^'^v;^ 1 ' (2 - 59) 

s ^ J 2\fWi m + \f¥ri Z2 - (2 - 60) 



Since the original quadrature components are uncorrelated, with variances (2.50), s\ and s 2 are 
also uncorrelated, with variances 



(Asl) = (Asl) = = (lAdf) . (2.61) 
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is indeed two vacuum units bigger than the added noise (2.11) of an 



The added noise, ^(g 2 + 1) 
ideal linear amplifier. 

In the high-gain limit, these two additional units of vacuum noise become irrelevant. On the 
other hand, it is easy to see how to convert this measurement-based model of linear amplification 
into an ideal linear amplifier. What needs to be done is to make the initial state \<f>) of modes b 
and c the vacuum state of d. Noticing that 



d = 5i(r) ® Sj(r)4=(-6 + c)S\(r) 5i(r) 
V2 



where the squeeze parameter is specified by 



•2r 



g-1 

5+1 



(2.62) 



(2.63) 



we realize that all we need to do is to use an initial state of the form (2.48), with r specified by 



this new value. The result is ideal linear amplification by a mechanism that is distinct from a 
parametric amplifier (see Fig. [4|: the primary mode has an input-output relation (2.58) that is the 
same as for a parametric amplifier; it is not hard, but tedious to check that the mode d, which is 
completely responsible for the added noise, does not evolve as does the single ancillary mode of a 
parametric amplifier. 

We can chase this new choice of initial state for the ancillary modes back through the above 
analysis to see what kind of measurement of x\ and X2 it corresponds to. The initial wave function 
for the ancillary modes becomes 



</3 2 ,Ti|0> 



1 



ire' 



-exp 



Pi + 7? 



2e 



2r 



(2.64) 



and this leads ultimately to the following Kraus operators for the simultaneous measurement of x\ 



and X2, replacing the coherent-state projectors (2.45): 



TT 



d 2 (3 



TT 



-(ff-l)[a-j9| s 



(2.65) 



This measurement gives up some sensitivity in determining the initial complex amplitude in return 
for maintaining enough coherence to introduce a bit less noise into the amplified output than does 



the Arthurs-Kelly measurement (2.45) 



5. Discussion 

Having surveyed various models for an ideal linear amplifier, we now turn to our main task, 
formulating a model of phase-preserving linear amplifiers and deriving the complete set of restric- 
tions on the noise that must be added in the amplification process. Our aim is to draw general 
conclusions, applicable to any phase-preserving linear amplifier. The standard input-output rela- 



tion (2.9), powerful though it is, is not sufficient for our purpose; the properties of the noise operator 



L are not sufficiently constrained, beyond the commutator that leads to the second-moment con- 



straint (2.11), to allow us to draw general conclusions about the full quantum statistics of L. We 
need a more precise characterization of the operation of a phase-preserving linear amplifier than 
the input-output relation. This we give in the next section. 

Before moving on, we note that Shi et al. |35j have considered two models of nonideal amplifiers, 
a laser amplifier with incomplete inversion and a cascade of alternating ideal amplifiers and ideal 
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attenuators, and found that the nonideal behavior of these amplifiers can be attributed to having an 
internal noise source that is not in its ground state. These findings provide additional motivation for 
our work and are consistent with our general conclusion that any phase-preserving linear amplifier 
is equivalent to a parametric amplifier with a physical initial state for the ancillary mode, an ideal 
amplifier arising uniquely in the case where the ancillary mode begins in the vacuum state. 



III. MATHEMATICAL CHARACTERIZATION OF A PHASE-PRESERVING LINEAR 

AMPLIFIER 

A phase-preserving linear amplifier takes an input signal to an output signal with the same 
phase, but with amplitude larger by a factor of the amplitude gain g. An essential feature of 
linearity is that the noise added by the amplifier is independent of the input signal. 

In this section, we capture this action mathematically in terms of two superoperators — linear 
maps on operators — that together characterize the operation of the amplifier. The first superop- 
erator accounts for the amplification by taking an input coherent state |a) to an output coherent 
state \ga): 

A(\a)(a\) = \ga)(ga\ . (3.1) 

The superoperator A amplifies without even the amplified input noise and so is clearly not physical 
by itself. The second superoperator includes the noise on the output signal by smearing out a 
phase-space distribution into a broader distribution: 

B = [ d 2 f3n^ 1 \(3)D(a,f3)QD\a,f3) . (3.2) 



Here marks the slot where the input to the superoperator goes. The real- valued function Il( _1 )(/3) 
is assumed to be normalized to unity on the phase plane; we call it the added-noise function 
(sometimes the smearing or spreading function). The added-noise function is independent of the 
input state, but it can and does depend on the gain g. We use a superscript (—1) on the added- 
noise function to indicate that II^ -1 ) has to do with an antinormal ordering, but the connection to 



antinormal ordering only becomes clear in Sec. IV A Other orderings for the spreading function 
will also arise as we proceed. 

The overall operation of the amplifier is given by acting first with A and then with B. This 
composition of the two is the amplifier map 

£(p)=B(A( P ))=BoA(p). (3.3) 

The natural operator-ordering perspective for the amplifier map becomes apparent when we deter- 
mine the output state for an input coherent state, 



£ 



[\a)(a\) = j d 2 prt- l \p)D(a^)\ ga )( ga \D\a^) = j d 2 (3 U^(p - ga)\ (3) 0\ . (3.4) 



For this input, the displaced spreading function, Il( -1 )(/3 — go), is the P function of the output 
state. Thus the perspective to have in mind is that of the P function in Fig. [2j A amplifies a 
coherent state without noise, and B, through the spreading function, accounts for all the noise at 
the output. 

The problem we are interested in can now be expressed as determining the restrictions on the 
added-noise function H"* -1 ) necessary and sufficient to ensure that the amplifier map £ can be 
implemented in a physical system. Mathematically, this is the requirement that £ be a completely 
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positive map, i.e., a (trace-preserving) quantum operation. We have already noted that A is 
not completely positive. Notice that B is a (trace-preserving) quantum operation if n( -1 )(/3) is 
nonnegative, but we do not need to assume this since it emerges from our main result |53j . 

Strictly speaking, a linear amplifier is phase preserving if and only if £ commutes with phase- 
space rotations. The amplifying map A does commute with rotations, but B commutes with 
rotations if and only if n( _1 )(/3) is independent of phase, i.e., depends only on \/3\. We do not 
need to assume that B has this property, however, to demonstrate our main result. The only 
phase-preserving property needed for our main result is built into A, i.e., that its raw amplification 
without noise is independent of phase. Thus we leave n^ _1 ^(/3) general for the present and make 
it independent of phase only when we consider examples of nonideal amplifiers in Sec. [V] and 
constraints on the moments of the added noise in Sec. IVI1 

An easily addressed point, which we use as an excuse to introduce a mathematical formulation 
we need, is that since the coherent-state projectors are not orthogonal, it is not obvious that A, as 



defined in Eq. (3.1), can be extended by linearity to all operators or, to put it differently, is even 
a linear map. To deal with this point, we actually define A by its action on the operator basis of 
displacement operators, which are <5-orthogonal: 

ti (D\a, a)D(a, /?)) = vr5 2 (a - /3) . (3.5) 

Here 5 2 (a) = 5(afj)d(aj) = 25{a\)5{a2) is the two-dimensional delta function on the phase plane. 
The definition of A becomes 

( 9 2 -l)|/3| 2 /2 9 2 

A(D(a, (3)) = - 2 D(a, (3/g) . (3.6) 



Linearity and the expression (12.46 ) for the coherent-state projectors as a Fourier transform of 



displacement operators can now be used to derive Eq. (3.1) as the action of the A on an input 
coherent state. 

These considerations suggest that it might also be useful to translate the action of B to the 
displacement-operator basis, 

B(D(a,f3)) = J d 2 aU ( -- 1) (a)D(a,a)D(a,p)D\a,a) = fl^* (p)D(a, (3) , (3.7) 

where 

fl(- 1 )(/3)= [ d 2 aU^{a)D(a^) = U { - 1) *(-^) (3.8) 



is the Fourier transform of the added-noise function. Normalization of the added-noise function is 
the statement that n^ _1 ^(0) = 1. Notice that these considerations allow us to write the amplifier 
map as£ = BoA = AoB', where B' is the same as B except that its added-noise function is 
nM)'^) = flt-V^/g) or, equivalently, U^'(a) = g 2 U^ l \ga). Since A o B' adds noise first 
and then does the amplification, this rescaling of the added-noise function is the map version of 
referring the noise to the input. 



Combining Eqs. (3.6) and (3.7) gives the action of £ in the displacement-operator basis: 



£{D(a,P)) = e^mV ^ ( - 1] y/9) D{ajm . (3 . 9) 



We can write the action ( |3.9[ ) more compactly in terms of antinormally ordered displacement 
operators, which absorb the Gaussian factors. A more general approach along these lines is to 
introduce the s-ordering of Cahill and Glauber |20t I21j. 

D^(a,f3) = e s|/3|2/2 £>(a,/3) , (3.10) 
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where s = gives symmetric ordering of products of a and a) , s = +1 gives normal ordering, 
-D( +1 )(a,/3) = e^ at e~^* a , and s = — 1 gives antinormal ordering, D^ 1 \a, (3) = e~ /3 * a e /3at . The 
s-ordered displacement operators satisfy the orthogonality relation 

tr(r>Wt( Q; a )Z)(-s)(a, 0)) = ^ 2 ( a - 0) . (3.11) 
In terms of s-ordering, Eq. ( |3,9| ) assumes the form 

£(£>«(a, /?)) = D( "f /g) flW(a, , (3.12) 

where we define an s-ordered version of IB" -1 ', 

= e ( S+ i) (g 2_ 1) ^|2 /2fi( _ 1) ^ ) = eS(fl 2_ 1)| ^ 2/2fi(0)(/3) ^ (3 l3) 

which satisfies nW*(/3) = nW(-jS) and nM(0) = 1. The corresponding (real-valued and normal- 
ized) s-ordered added-noise function is 

nM(a) = | ^n( s )(^)L>(/3,a) , (3.14) 

although we have no warrant that for s > — 1, this function is nonnegative or even exists. 
Our objective now is, first, to determine how the s-ordered characteristic function, 

$(*)(/?) = (D^(a,a))=ti{pD^(a,a)) , (3.15) 

transforms from an input state p to the state £(p) at the output of the linear amplifier and, 
second, to Fourier transform this result to find the corresponding input-output transformation of 
the s-ordered quasidistribution, 

wW(a) = / ^*«G8)I>GM) , (3.16) 



p 

where s = 1 gives the P function, s = gives the Wigner function, and s = — 1 gives the 
Q distribution. For this purpose, it is useful to translate Eq. (3.12) to the adjoint £*, which is 
defined by 

tr(Atf(B)) =tr([£*(A)]^B) (3.17) 

and which can be thought of as the Heisenberg-picture version of £. Using the 6- 
orthogonality (3.11), we have 

tr([£*(D^(a,a))]^D^ s \a,^)) =ti^D (s ^{a,a)£(D^{a,P))^ = U^ s > {a)nd 2 (/3 - go) . 

' ' (3.18) 
Comparing the leftmost and rightmost sides of this equality and again using Eq. (3.11) gives us 

£* (I>W (a, a)) = IT" 5 ) {a)D^ (a, got) . (3.19) 

We can now use the action ( |3.19| ) of £ * on displacement operators to determine the input-output 
transformation of the characteristic function: 

^:l(P) = tT{D^\a,P)£(p)) 

= tr([r(Z^)(a,-/3))]V) 
= &- s \p)ti{D^(a,gP)p) 

= u^(j3)^' ) { g p). (3.20) 
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The output s-ordered quasidistribution is a convolution of the s-ordered input quasidistribution, 
scaled by the gain, with the (— s)-ordered spreading function: 



[a 



W- 



W/9) 



r 



(3.21) 



The input-output transformation (3.21) is the generalization of the ball-and-stick depictions in 
Fig. [2] The P function of an input coherent state is a 5 function, so the output P function for 
this input, as noted above, is given directly by Il( _1 )(a), displaced to the position of the amplified 
expectation value of the complex amplitude. An ideal linear amplifier adds no nois e to the Q 
distribution, so H^ +l \a) = S 2 (a), which gives the input-output tranformation (2.27) for the Q 
distribution. The rescaled Wigner function of an input coherent state I7) is 



I e -2|/3-97| 2 /9 2 
7T 



9 2 vr g 2 /2 ' 

for an ideal linear amplifier, this rescaled Wigner function is convolved with 

, n . 1 e -2|a-/3| 2 /(9 2 -l) 

n(°V-/3) 

giving an output Wigner function 



tt (g 2 - l)/2 



1 e -|7-/?| 2 /(9 2 -l/2) 

Wout(a) = 2 77o 

vr g 2 - 1/2 



(3.22) 



(3.23) 



(3.24) 



which has the minimum output noise (2.13) permitted by quantum mechanics. 

We now want to go beyond these simple Gaussian considerations and to derive the general con- 
straints that complete positivity of £ places on the added- noise functions 11^. A straightforward 
approach invokes the Kraus representation theorem [29l [30] to conclude that if £ is a quantum op- 
eration, then there exists an ancilla E, with initial (pure) state \cp){4>\, and a joint unitary operator 
U such that 



£( P )=tT E {Up®\4>){<i>\tf) . 



(3.25) 



This is called an ancilla model for £ or a Stinespring extension 

It is useful to convert this ancilla model to the adjoint £*. Given any operators A and B on 
the primary mode, we have 

tx a ([£*{A)}^B) =tr a (A^£(B)) = tr ajE (rfUB ® {^(^tf) = tr a ((cf>\tf A^U\cf>)B) , (3.26) 

which implies, since B is arbitrary, that 

£*(A) = {<f>\U^AU\<j>) . (3.27) 

We can interpret the Kraus representation theorem as saying that £ is completely positive if and 
only if there exists a joint unitary and an ancilla state such that £* satisfies Eq. (3.27) for all 
operators A. 



In particular, using Eq. (3.19), we have that £ is completely positive if and only if 
n(-) (£)!>(-> (a, g(3) = £* {D^ (a,P)) = (<f>\ rfD& (a, 0)U\ 4>) . 



(3.28) 



This expression restricts the way U acts on joint states of the form (g> \(f>), where is an 
arbitrary pure state of the primary mode, and thus seems to promise a way to derive restrictions 
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on the added-noise function, but we have been unable to disentangle those restrictions from the 
freedom in choosing U and \(f>). Thus we drop this approach in favor of a different one (leaving 
open the question of whether the approach based on the Kraus representation theorem can be 
used to get to our main result). Instead of relying on the Kraus representation theorem to provide 
an ancilla model with a joint unitary and a physical ancilla state, we construct an explicit ancilla 
model using a particular joint unitary, the two-mode squeeze operator, and a "state" a of the single 
ancillary mode. The "state" a determines the added-noise function, and what we prove is that a 
must be a physical state, i.e., a valid density operator. 



IV. QUANTUM CONSTRAINTS ON ADDED-NOISE FUNCTIONS 
A. Two-mode squeezing model for any phase-preserving linear amplifier 

To develop this second approach, we begin with the s = version of the input-output transfor- 



mation (3.20) for the characteristic function. We define a unit-trace, Hermitian operator a, which 



we associate with an ancillary mode b, by 

d 2 /3 



a 



7T 



(4.1) 



By the completeness and orthogonality of the displacement operators, this is equivalent to 

m<J3) =tr b (D(b, sfcf^l /3>) = <^ 0) [V^IP*) , (4.2) 

where <J>(°) is the symmetrically ordered "characteristic function" of the "state" a. Quotes are used 
here because we have no warrant to assume that a is a valid density operator, i.e., has nonnegative 
eigenvalues. Indeed, the problem we are interested in, the restrictions on the added-noise function 
needed to ensure that S is completely positive, are translated to the corresponding restrictions 



on a. We show in Sec. IV B that a must be a valid state of the ancillary mode b. Formally, this 
means that a is a positive operator, having only nonnegative eigenvalues, a property denoted as 
a > 0. 



We can convert the characteristic function (4.2 ) to arbitrary ordering and then Fourier transform 



to relate the added-noise functions to the corresponding s-ordered quasidistributions of a: 

nW(/J) = *M(V7^1/3*) , (4.3) 

nW(a) _ wi-\-*/^p _ 



These definitions become useful when we introduce the two- mode squeeze operator (2.21) as a 



joint unitary operator for modes a and b; employing the input-output relation (2.22), we obtain 
tr a (D(a, P)ti b (Sp ® aS*j\ = tT a>b (S^D(a, P)Sp ® a) 



tT a (D{a,g/3)p)tv b [D(b, y/g*-\p*)a 

d> i ( n 0) ( 5 /3)n( )(/3) 

tr a (D(a,P)£(p)) ■ 



(4.5) 
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TABLE I: Linear-amplifier input-output transformations for characteristic functions, $o Ut and $ in , and 

(s) (s) 

quasidistributions, W ou j! and W- n . The input-output transformations are specified by added-noise functions, 
n( s ) or ]j( s ) , P, W, and Q denote the Glauber-Sudarshan P function, the Wigner W function, and the 
Husimi Q distribution. The expressions in square brackets specialize the added-noise-function formula 
immediately above to an ideal linear amplifier, i.e., to ancillary- mode state a - 



s arbitrary 



s = +1 



s = 



s = -1 



Characteristic-function transformation 

$il(^)=fK- s )(/3)^ ) (5/3) 

e -(l+ S )(g 2 -l)|/3| 2 /2" 

$i+ 1) (/3)=n ( - 1) (/3) < f i ( „ +1) (5/3) 
ft(- l >(0) = 4 _1) (^g^lf3*) 

e -(ff 2 -l)l/3| 2 

^ u ) t (^) = n(°)^)^° ) (^) 
n(°)(/3) = $( 0) (y^i>) 

e -(9 2 -l)l/3| 2 /2" 

C t 1) (/3) = n( +1 )(^)47 1) ( ff /3) 
n(+ 1 )(/3) = $i +1) ( v ^I / 3*) 
[i] 



Quasidistribution transformation 



W { - s) (-a* 



9 2 1 



s 2 -i 



1 e -2| Q | 2 /(l+ S )( 9 2 -l) 
* (l + S )( S 2-l)/2 



Po«t(a)= I d^-^a-P)^^ 



n H, (a) = ^WH 

.9 — 1 

1 e -l Q l 2 /(s 2 -i) 



7T ^-1 



w. n)= /^ n (o)(a-/3)^M 

^(- a y^T) 

2 //„2 



n(°)(a) = 



X e -2|a| 2 /(g 2 -l) 



7T (. 9 2 -l)/2 



, - / j2 fln (+l)/ a -,Qm[0/g) 
t(a) = / d /3IT^ ; (a - /?) ^ 



n< +1 )(a) = 



1 2_ 1 



[S 2 (a)] 



Since the displacement operators are a complete, <5-orthogonal set of operators, we get 

£(p) = tr b (Sp<g>aSt) , (4.6) 

just as though we had the ancilla model of Fig. [3jb), i.e., a single ancillary mode b that interacts 
with the amplifier mode a via a two-mode squeezing interaction. 

It is trivial that if a > 0, then £ is completely positive, so what we have to prove is that if £ is 
completely positive, then a > 0. Jiang, Piani, and Caves [55] have established the necessary and 
sufficient properties of a unitary operator U such that the complete positivity of a map 

£{p) = tr E (Up ® arf) (4.7) 

implies that the ancilla "state" a is a valid density operator. They show that a sufficient, but 
not necessary condition on U is that it be full rank, i.e., that the ancilla operators in its operator 
Schmidt decomposition span the space of operators on the ancilla. The two-mode squeeze operator 
is full rank, so we could rely on the results of [55] to assert our main result. Since the proof in [55] 
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assumes finite dimensions, however, we first prove, in Sec. IV B that S(r) is full rank for r / 
and then use this result to show that a > 0. 

Our main result thus is that any phase-preserving linear amplifier is equivalent to a parametric 
amplifier with a physical state a for the ancillary mode. An ideal linear amplifier is the case where 
a is the vacuum state. Stated in terms of the added-noise functions, our result shows that they are 
rescaled quasidistributions of the state a. In particular, for our formulation of a linear amplifier 
in terms of the maps A and B, which uses the P- function perspective of Fig. [2j the added-noise 
function Il( -1 ) is a rescaled Q distribution of the ancillary mode; the added-noise function Il( -1 ) is 
thus required to be everywhere nonnegative, although this is by no means sufficient to guarantee 
that the added-noise function is physical. If, instead, we use the Q-distribution perspective of 
Fig. [2J the added-noise function H"( +1 ) is a rescaled P function of the ancillary mode. Finally, if 
we use the symmetrically ordered moments that are usually used to discuss amplifier noise, the 
added-noise function n(°) is a rescaled Wigner function of the ancillary mode. These relations are 
summarized in Table |U 

We also note that for an ideal linear amplifier, the input-output transformation can be written as 
$W (£) = *W(gP) or, equivalents, as W^(a) = w£'\a/g)/g 2 , where s' = s- (1 + «)(1- 1/g 2 ). 
If, for example, s = 1, then s' = 1 — 2(1 — 1/g 2 ), implying that in the limit of high gain, the 
output P function looks like a rescaled input Q distribution. For arbitrary input states, as g — > 00, 
the amplifier noise wipes out any singularities or negativity associated with the input P function. 
Nonetheless, Nha, Milburn, and Carmichael [42] have shown that there are input states for which 
nonclassical features persist in the output for arbitrarily large gain. 

B. Proof of main result 

To prove our main result, we first establish that the squeeze operator is full rank; i.e., we show 
that given an operator Ob on the ancillary mode, tr^^O;,) = implies Ob = 0. We begin by 



writing the two- mode squeeze operator in a form similar to that in Eq. (2.24): 

S f r \ = g ~(ata+btb+l)/2 e -aTbtiy^ e ab^^l g -(a''a+btb+l)/2 _ ^ gN 

This gives us, for any coherent state \a) on mode a, 

(a\S(r)\a) = e -M 2 d^/ 9 ) g -(^ + i)/2 {a/ ^\ e -«W^ e ab^ la/V g )g -Vb+i)/2 

= e nH 9 -if/2 gg -^b + i)/2 D ^ _ ^ 2 -V V (btfe+1)/ " • ( 4 -9) 

Our premise, that tTb(SOb) = 0, implies that for all a, 

= tr b ((a\S\a)O b ) = e^-^/ 2 Hr b (d (b, -Ji—ta* ) g ^+l)/2 0bg -^b+i)/2 \ {4 W) 




Since the displacement operators are a complete, <5-orthogonal set of operators, we get that 
g -^b+i)/2 0bg -^b+i)/2 = o and) henc6) by the invertibility of g -^ b + l )l' 2 i that O b = 0. This 

establishes that S is full rank. 

Suppose now that we have a full-rank joint unitary U and a quantum operation defined as in 



Eq. (4.7). Diagonalize the ancilla initial "state" a, and decompose it into positive-eigenvalue and 



negative-eigenvalue parts, 

a = a + - cr_ , (4.11) 
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where 



°+ = J2 X j+\ e 3 + )( e j+\ > 



J+ 



°- = ^2 x j-\ e j-)( e j~ 



(4.12) 
(4.13) 



The eigenvectors \&j±) make up an ortho-normal basis. We assume that the Xj+s are strictly positive 
and allow the Xj~s to be positive or zero (zero so as to fill out the orthonormal basis with zero- 
eigenvalue eigenvectors). If we take the ancilla trace in a basis \fk), the quantum operation takes 
the form 



£{p) = M kd+ pMl j+ - M kJ .pM ) 

k,j+ k,j- 



k,j- ' 



where 



M kJ ± = y/^Z(f k \U\e j± ) = y/XjZtT b (U\e j± )(f k \) 



(4.14) 



(4.15) 



are "operation elements" that decompose £ into positive and negative parts. 

If any nonzero negative-part operation element, say M^ j_, does not lie in the operator subspace 
spanned by the positive-part operation elements, {M k j + }, then by projecting Mk,j- orthogonal 
to this operator subspace, we obtain a nonzero operator Nj^- such that 

El tr K*- M ^+)l 2 - El tr K^- M ^-)l 2 = - El tr K*- M ^-)l 2 

k,j+ k,j- k,j- 

< -\tv(N\ K _M J:K ^)\ 2 

< , (4.16) 

which implies that £ is not completely positive (for a completely positive £, which thus has a 
Kraus decomposition, there can be no operator like Nj t x-)- Thus, no matter what U is, complete 
positivity of £ requires that all the negative-part operation elements lie in the span of the positive- 
part operation elements. 

This means that for any K and J— , 



k,j+ 



■3 + ' 



(4.17) 



for some coefficients c? , J , . Using the definition (4.15), we can rewrite this expression as 



"k,j+ 



= tv b lu(^XjZ\ej-)(f K \-J2' 

V v kj+ 



c kf+ \Ai+l e i+)(A 



That U is full rank implies that 



^J\JZ\ej-){f K \ - %'/+ V x ~j+\ e j+)(fk\ = • 

k,j+ 

Applying this expression to \fi), we get 

= y/\JZ\ej-)5 K i - ^c5+"v / Vfl e i+> • 
3+ 



(4.18) 



(4.19) 



(4.20) 



24 



Since the vectors \ej±) are linearly independent, all the coefficients in this expression must be zero. 
In particular, for / = K, we get that Aj_ = 0. Since this holds for all J — , we can conclude that 
<j_ = 0. Thus, we reach the desired conclusion that for a full-rank U, S(r) being an example, a 
must be a valid density operator. 

That a must be a positive operator has been proven for the amplifier transformation and for more 
general phase-space transformations [56], using a technique that does not reveal the connection to 
a particular unitary operator, the two-mode squeeze operator in the case of a linear amplifier. 



V. EXAMPLES OF NONIDEAL LINEAR AMPLIFIERS 



We reiterate that the key assumption necessary for our proof is the phase-preserving amplifica- 
tion carried out by the map A. The proof does not require that the added-noise function be phase 
insensitive, i.e., that a be invariant under phase-space rotations, nor even that a have zero mean 
complex amplitude. Since phase-preserving amplifiers do have phase-insensitive noise, however, 
we assume that a is rotationally invariant for the examples considered in this section and for the 



discussion in Sec. VI of quantum limits on moments of the added noise. Rotational invariance 
implies that a is diagonal in the number basis 



\ n \n){n\ . 



(5.1) 



n=0 



In this section we consider both physical and unphysical a of this form, thus allowing us to compare 
physical and unphysical linear amplifiers. 



For rotationally invariant cr, the added-noise number ([2.18 ) is given by 

(\b\ 2 ) = (tfb) a + 



An 



1 



(5.2) 



The second- moment quantum limit, A\ > 1/2, becomes the constraint that a have a nonnegative 
mean number of quanta, (b%) a = {n) a > 0. 

The s-ordered characteristic function for a number state 1311 is 



4 s) (/3) 



2 / 2 (n\D(b,P)\n)=e( s - 1 ^ 2 / 2 L n (\tf), 



(5.3) 



where L n denotes the Laguerre polynomial of degree n. Fourier transforming gives the correspond- 
ing quasidistribution for a number state [3~T] . 



7T Z 



8+1 



7r(l — s) \S — 1 



-\0{- 



(5.4) 



We can plug these results into Eqs. (4.3) and (4.4) to find series rep rese ntations of the added-noise 



functions, IT S ) and n( s ), for the general rotationally invariant state (5.1 ). The series representation 
is particularly useful when a has only a few nonzero eigenvalues. 



In this subsection, we use the P-function perspective introduced in Sees. II B and III for which 
we need the Q function of a number state, 



1 



Q n {a) = -\{a\n)\ 



-H 2 \a\ 2n 



7T 



7T 



The added-noise function for the rotationally invariant state is 

Q 



n(^)( 



oo 

E 

n=0 



An" 



-a 



JL 
1 



1) 



o 



;-M 2 /(9 2 -l) °^ 



in 



(5.5) 



(5.6) 
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FIG. 5: Output P functions for five amplifiers. For all five, the input is a coherent state \f3), with (3 = 1, 
and the amplitude gain is g = 4; these correspond to the conditions in Figs. [T] and [2j The output P function 
is given by the added- noise function as in Eq. ( |5.8| . In (a) the ancillary-mode state a is the vacuum state, 
which gives an ideal linear amplifier. The inset shows the P function of the input coherent state as a tiny 
black dot and the output P function, P ou t(aO = e~' Q ' ^ 9 /n(g 2 — 1), as colored contours; this inset is the 
full contour plot of the P-function perspective that is given as a stick figure in Fig. [2j The main plot in (a) 
displays the same output P function in a thre e-dimensional plot. In (b)-(e), the ancillary-mode "state" is 
a = (| - A)}0><0j +A|1)(1| + i]2><2j [Eq. fcti ]: (b) A = 0.5, (c) A = 0, (d) A = -0.5, and (e) A = -1.0; the 
output P functions, given by Eqs. ( |5.8| ) and ( |5.9[ ), are displayed as three-dimensional plots. Plots (b) and (c) 
are physical amplifiers, whereas (d) and (e) are unphysical, this despite the similarity, at least by eye, of (d) 
to the Gaussian of the ideal linear amplifier in (a). All four values of A are such that the second-moment 
quantum limit is satisfied, which requires A > — 1, and the output P function is everywhere nonnegative, 
which requires -1.37 = -(1 + V3)/2 < A < 1/2. 



To illustrate the possibilities for the added noise, we specialize now to a one-parameter family 
of states, which have support only on the first three number states, 



a 



A |0)(0| + Ai|l)(l| + A 2 |2)(2| = (i - A)|0)(0| + A|l)(l| + §|2)(2| . (5.7) 



For this a to be physical, the three eigenvalues must be nonnegative, which restricts the parameter 
A to the range < A < i; here we will also be considering values outside this range, which give us 
unphysical a and, hence, unphysical amplifiers. The mean number of quanta, (b^b) a = Ai + 2A 2 = 
A + 1, gives a second- moment constraint, A > — 1; that this second- moment quantum limit allows 
unphysical values of A indicates that quantum constraints on higher moments are important. 
In our examples, we assume that the input to the amplifier is a coherent state In this situ- 



ation, the output P function is given by the added- noise function (5.6), displaced to the amplified 
mean complex amplitude, i.e., 

P out (a) = n(" 1 )(a- 5 /3) . (5.8) 
We can ask for the range of values of A for which the added-noise function, 

|2 /f„2_ 



e -MV(9 2 -l) / U|2 I | a |4 
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FIG. 6: Output P function (5.8), which is a displacement of the added-noise function (5.9), for amplifiers 
with ancillary-mode "state" a of Eq. (5.7). The parameter A varies from 0.5 to —1.5 in steps of size 0.025. 
As in Fig. |4j the input state is a coherent state \(3) with (3 = 1, and the amplitude gain is g = 4. The P 
function is plotted along the real a axis, which is sufficient because of the rotational symmetry of the noise. 
The plotted P functions are not normalized to unity; instead they are scaled to have maximum value equal 
to 1. The amplifiers in the range < A < 0.5, all of which are physical, are plotted in green and blue; 
A = 0, the boundary between physical and unphysical, is plotted as a black dashed line. When A = 0.5, 
there is no vacuum contribution to a, so the P function has a deep hole in the middle. As A decreases from 
0.5, the vacuum contribution increases, and the hole disappears, to be replaced by a distribution with broad 
shoulders, still faintly evident at A = 0. The transition from distributions that have a minimum in the 
middle (blue) to those with a maximum in the middle (green) occurs at A = 0.22. The values in the range 
— 1 < A < 0, in which the second- moment constraint is satisfied (A\ > 1/2), are plotted in purple; A = — 1 
is plotted as a black dashed line. The values A < —1, plotted in red, violate the second-moment constraint, 
and for A < — (1 + \/3)/2 = —1.37, the output P function takes on negative values. The Gaussian output P 
function of an ideal linear amplifier is plotted as a solid line; it reveals the nonGaussian character of all the 
added-noise functions associated with the ancillary-mode "state" a of Eq. (5.7). 



is everywhere nonnegative. We need Ao > to keep the added-noise function nonnegative at a = 0, 
and we need A2 > to keep it nonnegative when |ev| 2 is large (if A2 = 0, we also need Ai > 0). 
When A2 > 0, the polynomial Ao + X\x + ^A2^ 2 , which multiplies the Gaussian in the added-noise 
function, with x = \a\ 2 /(g 2 — 1), has a minimum at x = — A1/A2, where it takes on the value 
Ao — A 2 /2A2- We don't care about this minimum when Ai > 0, since the minimum then occurs 
at a negative value of x, but when Ai < 0, ensuring that the added-noise function is nonnegative 
requires that Ao — A 2 /2A2 > 0, i.e., 2A0A2 > A 2 . The requirements for nonnegativity of the added- 
noise function are thus (i) A > 0, (ii) A 2 > 0, and (iii) 2A A 2 > A? if Ai < 0. Translated to the 
single parameter A, the requirements for the nonnegativity of the added-noise function reduce to 
-(1 + y/3)/2 < A < 1/2. 

Figure [5] d isplays the output P function (5.8), which is a displacement of the added-noise 
function (5.9), for four values of A, two physical and two unphysical, using the ancillary-mode 
state (5.7). For comparison, the figure also shows the output P function for an ideal linear amplifier 
(Ao = 1, o = |0)(0|) with the same input state and the same gain. The lessons to be drawn from 
Fig.[5]are, first, that eyeballing the added-noise function is not a reliable way to assess its physicality 
and, second, that the second-moment quantum limit is not sufficient to discriminate physical from 
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unphysical amplifiers, since all four values of A satisfy the second- moment constraint. Figure [6] 
provides a more detailed look at the output P function for values of A between 0.5 and —1.5. 
The Gaussian output of an ideal linear amplifier is shown for comparison, making it easy to see 
the nonGaussian character of the output noise. Even so, it is not easy to judge whether these 
added-noise functions correspond to physical amplifiers just by looking at the plots, except for 
those values, A < —(1 + v3)/2 ; for which the output P function goes negative. 



VI. QUANTUM LIMITS ON ADDED-NOISE MOMENTS 

In this section we turn our principal result into quantum constraints on the moments of the 
added noise. We assume that all the noise, both the noise carried by the input signal and the added 



noise, is phase insensitive; for the added noise, this implies that a has the form (5.1). The only 
nonzero moments of Aa = a — (a), which we call noise moments, are those for which the number 
of creation operators matches the number of annihilation operators. 

We characterize the noise in terms of symmetrically ordered noise moments. For this purpose we 
introduce the notation | a \ 2k for the symmetric product of k annihilation operators and k creation 
operators; formally, we can write 



, 2fc _ d 2k D(a, a 



da k d(—a 



*\k 



,n 2 



(2fc)! 



sum of the (2k)\/(k\) 2 ways of ordering 
a product of k annihilation operators | • (6-1) 
and k creation operators 



In terms of this notation, the nonvanishing input and output noise moments are written as (|Aa| 2fc ). 
Details of our moment manipulations are relegated to an Appendix. 
The input-output relation for the nonvanishing noise moments is 

(\Aa out n = £ ( , )V' m V - ir<|Aa in | 2 (<— V m , (6.2) 

m=0 ^ ' V 

where 



M 2fc L, (6-3) 



the (2fc)th symmetric moment of b and , we call the kth added-noise number. Equation (6.2) 
expresses how the the input noise, given by the moments (| Aai n | 2 ( fc-m )), combines with the added 
noise, given by the added noise numbers A m , to produce an output noise moment. The last (m = k) 



term in the sum (6.2 ) comes only from the added-noise number A^; it characterizes the noise added 
at the (2fc)th moment. If the signal noise is known a priori, then the added-noise numbers can 
be obtained from measurements of successive output noise moments; such a procedure has been 
implemented, using dual input and output ports, in 



Using Eq. (A.6), the added-noise numbers can be written in terms of normally ordered 



moments, or factorial moments, ((tf) m b m ) a = {(tfb)(tfb - 1) • • • (tfb - m + 1)) CT , 

A " = I £ uy J -l^mc ■ (6.4) 
2 K ^-^ muk — m)\ m\ 



Since all the terms in this sum are nonnegative, A^ is bounded below by the state-independent 
m = term, which is equal to k\/2 k . This bound is achieved if and only if all the higher terms in 
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the sum vanish, and that occurs if and only if a is the vacuum state. Thus we have the following 
quantum limits on the added- noise numbers: 



^><0||6| 2fc |0) 



For comparison, a thermal state, 



1 



a 



n 



1 + n \ 1 + n 



k\ 

fetfc 



with mean number of quanta n, has factorial moments ((b^) k b k ) 
A k = k\(n + \f. 



(6.5) 



(6.6) 



k\ n and added- noise numbers 



The quantum limits (6.5) are not very useful, except for the familiar second- moment constraint 



Ai > 1/2. One way to see this is to return to the family of "states" of Eq. (5.7). For any rotationally 
invariant a, as in Eq. (5.1 ), for which A n = for n > N, the factorial moments ((b^) k b k ) a vanish for 
k > N. For the states of Eq. (5.7), all the factorial moments vanish except (b^b) a = A1 + 2A2 = A + l 
and ((6 t ) 2 6 2 ) (T = 2A 2 = 1, and these give added-noise numbers A k = (A;!/2 fc )(l + 2/c[Ai + (/c + l)A 2 ]). 
The quantum limit (JOj) is satisfied if Ai > —(k + 1)A 2 [A > —(k + l)/2]. Thus the "states" with 
A > —1 satisfy the quantum limits (6.5) for all k; in particular, the two unphysical "states" depicted 



in Fig. [5] satisfy all these quantum limits. 

To do better, we need conditions on the noise that are necessary and sufficient to guarantee 
that (j is a valid density operator. Since we specialize to phase-insensitive added noise, for which 
a is diagonal in the number basis, we are dealing with a classical probability distribution, defined 
by the eigenvalues A n . Thus the appropriate conditions can be obtained from the solution of the 
classical moment problem [57-59J: what sequences of moments are consistent with a (necessarily 
nonnegative) probability distribution? The answer to this question depends on the domain on which 
the probability distribution is defined: the Hamburger moment problem deals with a distribution 
defined on the entire real line, i.e., —00 < x < 00; the Hausdorff moment problem deals with a 
distribution defined on a finite interval, which can be taken to be < x < 1; neatly sandwiched 
between these two problems, the Stieltjes moment problem deals with the domain < x < 00 and 
thus provides the meat to feed into the maw of our analysis. 

The moments in our situation are moments of the number operator: 



M l = {{b%) l ) a = (n l ) a 

The Stieltjes problem is stated in terms of a sequence Mi, I 
two sequences of matrices, 



1,2, 



(6.7) 

from which one constructs 
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M 3 M 4 • • • 


M k+2 
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•• M 2k+ J 



(6.8) 



(6.9) 
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k = 0,1, 2,.... The sequence Mi consists of moments of a nonnegative probability distribution if 
and only if (i) det Q { ° ] > and det > for k = 0, 1, 2, . . . or (ii) det Q { ° ] > and det qJ. 1} > 

for k = 0, 1, . . . , K and det = and det = for k > K . Case (i) leads to a distribution 
with infinite support, and case (ii) to a distribution with finite support. These then are the 
amplifier quantum limits expressed in terms of number moments. In the following, we give these 
quantum limits as in case (i), i.e., with strict inequalities, but discuss the consequences of case (ii) 
for equalities in the quantum limits. We take no account of the fact that we are concerned with 
distributions concentrated on the nonnegative integers, whereas the Stieltjes problem deals with 
distributions on the continuous domain < x < oo. 

The first four nontrivial quantum limits imposed by the solution to the Stieltjes problem are 
the following: 

< Mi , (6.10a) 
< M 2 - Ml , (6.10b) 
0<MiM 3 -Mf, (6.10c) 
< M 4 (M 2 - Ml) - M 3 (M 3 - MM) + M 2 (M 1 M 3 - M 2 2 ) . (6.10d) 

The consequence of case (ii) is that there can be equalities in this list, but once one encounters an 
equality, all subsequent constraints in the list must also be equalities. 



The first constraint (6.10a) is simply that the mean number of quanta is positive. Allowing for 
equality in this constraint gives the usual second- moment quantum limit. Equality implies that 
a is the vacuum state (thus an ideal linear amplifier); since all higher moments also vanish, all 



the constraints become equalities. The second constraint (6.10b) requires that the variance of the 
number of quanta be positive. The consequence of equality in this constraint is that the variance 
is zero, which implies that a is a number eigenstate and thus that all the constraints except the 
first are equalities. Notice that the first three constraints imply < Mi(M 3 — MiM 2 ), which since 
Mi > 0, implies < M 3 - MiM 2 . 



We can apply the number- moment quantum limits (6.10) to the one-parameter family of 
"states" (5.7), for which the number moments are M; = Ai + 2'A 2 = A + 2 l ~ 1 . The first four 



quantum limits reduce to A > -1; A 2 + A - 1 < 0, i.e., -(l + \/5)/2 < A < (— l + VE)/2; A > 0; and 
A(| — A) > 0. The first two quantum limits do not rule out all unphysical states, the third rules 
out unphysical states with negative values of A, and the fourth rules out all unphysical states. One 
cannot achieve equality in the first two constraints, since doing so is unphysical, but the third and 
fourth constraints can achieve equality. 

To be able to use the number- moment quantum limits, we need the relations between number 



moments and the added-noise numbers. Using Eq. (A. 13), the added-noise numbers can be written 
in terms of number moments as 

h h 

k\ -A k\ 



1=0 



m\{k — m)\ m 



where denotes a (signed) Stirling number of the first kind [60], EI]. The first four cases of 



Eq. (6.11) are the following: 



Ai = Mi + i , (6.12a) 

,42 = M 2 + Mi + i, (6.12b) 

A 3 = M 3 + |M 2 + 2Mi + f , (6.12c) 

J 4 4 = M 4 + 2M 3 + 5M2 + 4Mi + | . (6.12d) 
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The constant term in these expressions is the quantum limit (6.5). More useful is to write the 



number moments in terms of added-noise numbers. The number moments can be written in terms 



of factorial moments using Eq. (A. 19) and in terms of the added-noise numbers using Eq. (A.20): 

I 



Mi 



k=0 



I 

£ 

m=0 



A r . 
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k=m 



k—m 



k\ 



k\ 



Jk) 



m\ m\(k — m)\ 



(6.13) 



Here S, . , defined in Eq. (A.16), denotes a Stirling number of the second kind \60\ 161] . The first 



four cases of Eq. (6.13) are the following: 



Mi 

M 2 
M 3 
M 4 



A X 

A 2 
A 3 
A 4 



i 

2 ' 

A x , 

\A 2 - \A X + \ 



2A, 



2--i ' 4 

2A 2 + 2Ai . 



(6.14a) 
(6.14b) 
(6.14c) 
(6.14d) 



Plugging these expressions into the number-moment quantum limits (6.10) gives the first four 
quantum limits in terms of the added-noise numbers: 



Ax > 



2 ' 



A 2 > A\ + 



1 



A 3> 7, [3A 2 +Ax- 
A A > 2{A 3 + A 2 -Ai 



+ 



(A 2 -A^ 



A, 



+ 



(A 2 - Axf + ±(4A 3 - 6A 2 - 2Ax + l)[8Ai(Ai - A 2 ) + AA 3 - 2A 2 - 6Ai + 1] 



Ao 



A\ 



(6.15a) 
(6.15b) 
(6.15c) 

(6.15d) 



The complexity of the last of these expressions suggests that the best way to deal with quantum 
limits on higher moments of the added noise is to translate measured added-noise numbers into 



number moments, as in Eqs. (6.10). 



effective number moments using Eq. (6.13) and to use the quantum limits expressed in terms of 



VII. CONCLUSION 

Amplification, by translating from the real world of quantum physics to the mundane world of 
everyday experience, is a principal means by which we gain access to the quantum world. Phase- 
preserving amplification transforms signals too weak to be perceived into much larger signals that 
we can lay our grubby, classical hands on. As the signal transitions to the classical world, however, 
quantum mechanics extracts its due: any phase-preserving linear amplifier must add noise, which 
is equivalent to half a quantum at the input in the limit of high gain. 

In this paper we consider the full set of quantum constraints on the operation of a single-mode 
phase-preserving linear amplifier, going well beyond the usual emphasis on second moments and 
Gaussian noise. Our main result is that any phase-preserving linear amplifier is equivalent to a 
parametric amplifier with a single ancillary mode that begins in a physical state a. The noise 
properties of the amplifier, even should it bear no resemblance to a parametric amplifier, are 
encoded in the effective state a. In particular, the noise the amplifier adds to a signal, as encoded 
in symmetrically ordered moments, is described completely by the Wigner function of a. Using this 
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general characterization of linear amplification, we consider how the phase-space quasidistributions 
of the noise input to the amplifier and the noise added by amplification combine to produce the 
noise at the output of the amplifier, and we derive quantum limits on all moments of the added 
noise. 

Despite the length of this paper, there is work still to be done. Perhaps the most important 
extension of our work will be to amplification of continuous-time signals. Such signals are best 
dealt with in the Fourier domain, where we can think of a phase-preserving linear amplifier as one 
that amplifies a continuum of frequency modes with a frequency-dependent gain. We expect our 
main result to generalize in the obvious way: at each frequency, the amplifier will be equivalent to 
a parametric amplifier, with an ancillary mode that provides the frequency-dependent gain at that 
frequency; the joint state of all the ancillary modes will have to be physical, but the ancillary modes 
will not have to be independent, even in the case of time-stationary noise. The second moments 
of the added noise will be expressed in terms of a spectral density of added noise, which will obey 
the usual quantum limit |6j|33]. For Gaussian noise, the added-noise spectral density will be the 
entire story, but for a nonGaussian amplifier, there will not only be the possibility of nonGaussian 
ancillary-mode states, as in a single-mode amplifier, but also the possibility of entanglement among 
the ancillary modes at different frequencies. 

A second extension involves how best to characterize the performance of a phase-preserving 
linear amplifier. We derive in Sec. VI the quantum limits on the measured moments of the added 
noise. These limits are both cumbersomely complicated and not really the point. The best way 
to characterize the performance of a linear amplifier would be to translate the measured noise 
into an effective ancillary- mode state a, i.e., into estimates of the eigenvalues \ n of a. A nearly 
quantum-limited amplifier, for example, will have Ao close to 1. Thus what one would like to do is 
to perform indirect tomography [55], in which one uses measurements on a system, in this case the 
amplified modes, to reconstruct the state of (perhaps imaginary) ancillas, in this case the ancillary 
modes of our parametric-amplifier model. In the amplifier context, this sort of tomography is a 
species of optical homodyne or heterodyne tomography [62, 63J, since one uses the statistics of 
linear measurements at the output of the amplifier, like homodyne or heterodyne measurements, 
to reconstruct a quantum state, in this case the state of the imaginary ancillary modes. This sort 
of tomography is a tricky business, fraught with instabilities. We defer consideration of it to future 
work, not just because it's tricky, although that is a problem, but also because the job really should 
be done in the context of continuous-time, multi-frequency linear amplification. 
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Appendix: Symmetrically and normally ordered products and number powers 



In this Appendix, we give the relations among ordered products and powers of the number 
operator that are used in Sec. VI For details on the hypergeometric function 2^1, see [60J, and for 
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details on the Stirling numbers, see [60| I6T]. 

The input-output relation for a phase-preserving linear amplifier, expressed in terms of sym- 



metrically ordered characteristic functions (s = 0), combines Eqs. (3.20) and (4.3) into 

*£l(# = *£W£ 0) (>/^/r). 



(A.l) 



If we assume that a has no mean field, as is true for the rotationally invariant a of Eq. (5.1), 
the expectation value of the complex amplitude of the primary mode transforms as in Eq. (2.7). 
Factoring out the input and output expectation values from the characteristic functions gives new 
characteristic functions, 



W) = e-M^'+P^^tP) = (D(Aa in ,f3)) , 
= e^ {a ° ut >* + ^° ut >$2(/3) = (£>(Aa out ,/3)> , 



(A.2) 
(A.3) 



which generate symmetrically ordered noise moments, i.e., moments of Aa = a — (a). In terms of 
these new characteristic functions, the input-output relation ( A.l[ ) becomes 



^°u ) t (/?) = <° ) (5/3)4 0) (V^/3*). 



,(o) 



(A.4) 



If we further assume that all the noise is phase insensitive, i.e., that the characteristic functions in 
Eq. (A.4) depend only on |/3|, then the only nonzero moments are those with an equal number of 
creation and annihilation operators. The input-output relation for these noise moments is 
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(A.5) 



In the second line, all the other possible derivatives vanish as a consequence of phase insensitivity, 
i.e., because the characteristic functions depend only on the absolute value of their arguments. The 
last term in the sum (m = k) comes only from the added noise and characterizes the noise added 
at the (2&)th moment, so we define it, in Eq. (6.3), to be the fcth added-noise number. 

To derive quantum limits on the added-noise numbers, we need to relate them to moments of 
the number operator b^b. We do this in two steps. The symmetrically ordered product (6.1) is 
related to normally ordered products by 
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(A.6) 



Notice that (n\(tf) k b k \n) = n{n- 1) • ■ ■ (n- 
k-l) = (x + k-l)\/(x-l)\, k = 0,1,2,.. 



- k + 1) = (— l) fc (— n)k =, where (x)k = x(x + 1) • • • (x + 
, denotes the Pochhammer symbol (or rising factorial). 
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The falling factorial can be written in terms of the Pochhammer symbol as (— l) fc (— x)k = x ( x ~ 
1) • • • (x — k + 1) = x\/(x — k)\ . In terms of this notation, the normally ordered product, 



(&t)*6* = 6+6 ( b t 6 _i)... ( b t & _ fc + i) = (_!)fc(_ 6 t & ) fc > 



(A.7) 



is the falling factorial. 

The second step is to write the normally ordered products in terms of powers of the number 
operator. One way to do this is to iterate the recursion relation, 



(&t&) Jfe = (&t&-fc + l)(&t&) fc _ 1 , 
to generate the required relations, the first four of which are 

6 f 6 = 6 f & , 

(&t)2&2 = (6 f 6 )2 _ b ] h ; 

(6t)3 & 3 = ( 5 t 6 )3 _ 3 ( 6 t & )2 + 2h ] h ; 

(ftt)4 6 4 = (6 f 6 )4 _ 6 ( 6 t 6 )3 + n(6 t 6 )2 _ Q b \ b _ 

Generally, we can use the expansion of the falling factorial as a polynomial in x, 
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(-i) fc (-*) fc = ^V, 

1=0 



(A.8) 



(A.9a) 
(A.9b) 
(A.9c) 
(A.9d) 



(A.10) 



where the coefficients Sl \ I < k = 0, 1, . . . , are the (signed) Stirling numbers of the first kind. The 

Stirling numbers satisfy = 8ko, which makes (x)o = 1. Equation (A. 10) converts normally 
ordered products to powers of the number operator: 



(6t) fc 6 fc = (_i)*(_6t & ) fc = £ 5 (0( 6 t b )J . 



(A.ll) 



Using the Pochhammer symbol, we can rewrite Eq. (A. 6) as 
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(A.12) 

where a-^i denotes the hypergeometric function. Plugging Eq. (A.ll) into Eq. (A.6) gives us a 
closed- form expression for symmetric products in terms of powers of the number operator: 
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The last term in each expression is equal to k\/2 k . 



We can invert Eq. (A. 13) by following the same steps in the opposite direction. The normally 



ordered products are related to symmetric products by 
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are the matrix inverse of the Stirling numbers of the first kind, i.e., 
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This can be used to invert Eq. (A. 10) 



k=0 



and, hence, to invert the corresponding operator relation (A. 11): 
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Now, plugging Eq. (A. 15) into Eq. (A. 19) gives us 
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